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ABSTRACT 

We present a pulsational stability analysis of hot post-AGB H-deficient pre-white dwarf stars with 
active He-burning shells. The stellar models employed are state-of-the-art equilibrium structures rep- 
resentative of PG1159 stars derived from the complete evolution of the progenitor stars, through 
the thermally pulsing AGB phase and born-again episode. On the basis of fully nonadiabatic pul- 
sation computations, we confirmed theoretical evidence for the existence of a separate PG1159 in- 
stability strip in the logT c ff — \ogg diagram characterized by short-period g-modes excited by the 
e-mechanism. This instability strip partially overlaps the already known GW Vir instability strip 
of intermediate/long period <?-modes destabilized by the classical K-mechanism acting on the partial 
ionization of C and/or O in the envelope of PG1159 stars. We found that PG1159 stars characterized 
by thick He-rich envelopes and located inside this overlapping region could exhibit both short and 
intermediate/long periods simultaneously. 

As a natural application of our results, we study the particular case of VV 47, a pulsating planetary 
nebula nucleus (PG1159-type) particularly interesting because it has been reported to exhibit a rich 
and complex pulsation spectrum including a series of unusually short pulsation periods. We found that 
the long periods exhibited by VV 47 can be readily explained by the classical K-mechanism, while the 
observed short-period branch below w 300 s could correspond to modes triggered by the He-burning 
shell through the e-mechanism, although more observational work is needed to confirm the reality 
of these short-period modes. Were the existence of short-period g-modes in this star convincingly 
confirmed by future observations, VV 47 could be the first known pulsating star in which both the 
K-mechanism and the e-mechanism of mode driving are simultaneously operating. 
Subject headings: stars — pulsations — stars: individual: VV 47 — stars: interiors — stars: evolution 
— stars: white dwarfs 



1. INTRODUCTION 

At present, the study of stellar pulsations constitutes 
one of the most fundamental pillars on which the building 
of stellar astrophysics rests on. Although the theory of 
stellar pulsations was initially elaborated to explain the 
existence of classical variable stars such as Cepheids and 
RR Lyrae, in the last few decades the study of pulsating 
stars has been extended to many other different kinds of 
stars that were regarded as constant stars before (e.g., 
Unno et al. 1989; Gautschy & Saio 1995). Nowadays, 
new classes of pulsating stars across the HR diagram 
are being routinely uncovered from ground-based obser- 
vations as well as space missions (e.g., CoRoT, MOST; 
see Aerts et al. 2008). The study of stellar pulsations 
through the approach of asteroseismology constitutes a 
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powerful tool to probe the internal structure and evolu- 
tion of stars. 

Most of the pulsations exhibited by pulsating stars are 
self-excited through the classical K-mechanism operating 
in a partial ionization zone near the surface of stars (Cox 
1980; Unno et al. 1989). As a matter of fact, this mecha- 
nism is responsible for pulsations of the stars in the clas- 
sical instability strip due to partial ionization of H and 
Hel and/or Hell. In the driving zone, the opacity pertur- 
bation increases outward so that radiative luminosity is 
blocked in the compression phase of pulsation. The re- 
gion gains thermal energy in the compression phase and 
it loses thermal energy in the expansion phase. 

A less common — and consequently less explored — pul- 
sation driving mechanism in stars is the e-mechanism. 
This mechanism is due to vibrational instability induced 
by thermonuclear reactions. In this case, the driving 
is due to the strong dependence of nuclear burning on 
temperature. During maximum compression, the tem- 
perature and thus the nuclear energy production rates 
are higher than at equilibrium. So, in the layers where 
nuclear reactions take place, thermal energy is gained at 
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compression while the opposite happens during the ex- 
pansion phase (Unno et al. 1989; Gautschy & Saio 1995). 
An excellent historical account of studies on vibrational 
destabilization of stars by the e-mechanism, can be found 
in Kawaler (1988) — we refer the interested reader to that 
paper for details. 

In this paper, we explore the e-mechanism in connec- 
tion with pulsating PG1159 stars. These stars, also 
called GW Vir or DOV, are very hot H-deficient post- 
Asymptotic Giant Branch (AGB) stars with surface lay- 
ers rich in He, C and O (Werner & Herwig 2006) that 
exhibit multiperiodic luminosity variations with periods 
ranging from 300 to 6000 s, attributable to non-radial 
g-modes driven by the K-mechanism acting on the re- 
gion of partial ionization of C and O (Starrfield et al. 
1983, 1984, 1985; Gautschy 1997; Quirion et al. 2004; 
Gautschy et al. 2005; Corsico et al. 2006; Quirion et 
al. 2007). Some pulsating PG1159 stars are still embed- 
ded in a nebula and arc called Planetary Nebula Nuclei 
Variable (PNNV) stars (see Winget & Kepler 2008 and 
Fontaine & Brassard 2008 for recent reviews). 

Evolutionary models of PG1159 stars with thick He- 
rich envelopes located at the upper left portion of the 
HR diagram are characterized by the presence of vigor- 
ous He-burning shells. The first attempt to study the ef- 
fect of the e-mechanism induced by a He-burning shell in 
H-dcficient pre-white dwarf stars was the seminal work 
by Kawaler et al. (1986). These authors found some 
g-modes excited through this mechanism with periods 
in the range 70 to 200 s. Observationally, however, no 
signature of these short pulsation periods was found in 
the surveys of planetary nebula nuclei conducted at that 
time (Grauer et al. 1987; Hine & Nather 1987). Later 
on, stability analysis on simplified PG1159 models by 
Saio (1996) and Gautschy (1997) also predicted unstable 
g-modes driven by the e-mechanism with periods in the 
range 110 — 150 s. 

The interest in the e-mechanism in the context of H- 
deficient post-AGB stars has recently been renewed by 
the discovery of luminosity variations in the PNNV star 
VV 47 (T cff = 130 000 ±5000 K, log<? = 7 ±0.5, C/He= 
1.5 and 0/He= 0.4; Werner & Herwig 2006) by Gonzalez 
Perez et al. (2006). The most outstanding feature of 
VV 47 is that its period spectrum appears to include a 
series of unusually short pulsation periods (~ 130 — 300 
s), the shortest periods ever detected in a pulsator of its 
class. Gonzalez Perez et al. (2006) (see also Solheim et 
al. 2008) speculate that these rapid oscillations could be 
excited by the e-mechanism. 

In this work, we largely extend the pioneering work by 
Kawaler et al. (1986), Saio (1996) and Gautschy (1997) 
by performing fully nonadiabatic pulsation computations 
on realistic PG1159 models extracted from full evolution- 
ary sequences with a wide range of stellar masses and 
effective temperatures. In particular, we gather strong 
evidence for the existence of a new short-period g-mode 
instability strip of pulsating PG1159 stars due to the e- 
mechanism. In addition, we examine the possibility that 
the short-period g-modes of VV 47 could be excited by 
this mechanism. The paper is organized as follow: in the 
next section we briefly describe the input physics of the 
PG1159 evolutionary sequences analyzed and the nonadi- 
abatic treatment of the pulsations. In Sect. [3]we describe 
the stability analysis. In Sect. 13.31 we present the appli- 




FlG. 1. — The PG1159 evolutionary tracks of Althaus et 
al. (2005), Miller Bertolami & Althaus (2006) and Corsico et 
al. (2006), with stellar masses of (from right to left): M* = 
0.515,0.530,0.542, 0.565, 0.589, 0.609, 0.664, O.741M . Also shown 
is the location of known PG1159 stars. The error bars for VV 47 
are displayed. 

cation to the star VV 47, and in Sect 13.41 we discuss the 
case of the prototypical DOV star PG 1159—035 in the 
context of our theoretical findings. Finally, in Sect. |4]we 
summarize our main results and make some concluding 
remarks. 

2. EVOLUTIONARY/PULSATIONAL MODELLING OF 
PG1159 STARS 

The PG1159 equilibrium models on which the present 
investigation rests on were extracted from the evolution- 
ary calculations presented by Althaus et al. (2005), 
Miller Bertolami & Althaus (2006) and Corsico et al. 
(2006), who computed the complete evolution of model 
star sequences with initial masses on the ZAMS in the 
range 1 — 3.75M Q . The evolutionary tracks for the H- 
deficient pre-white dwarf remnants are displayed in Fig. 
[TJ All of the post-AGB evolutionary sequences, com- 
puted with LPC0DE (Althaus et al. 2005), were followed 
through the very late thermal pulse (VLTP) and the 
resulting born-again episode that gives rise to the H- 
deficient, He-, C-, and O-rich composition characteristic 
of PG1159 stars. For details about the input physics and 
evolutionary code used, and the numerical simulations 
performed to obtain the PG1159 evolutionary sequences 
employed here, we refer the interested reader to the works 
mentioned above. One distinctive feature, and crucial 
for this study, that is common to all of our sequences, 
is that the PG1159 models are characterized by He-rich 
envelopes thick enough as to sustain active He-burning 
shell sources during the evolutionary stages of interest. 
This is at variance with the non-standard PG1159 mod- 
els employed in Althaus et al. (2008) to explain the II 
values in PG 1159—035, which are characterized by thin 
He-rich envelopes and so they are not able to sustain an 
active He-burning shell. 

The pulsational stability analysis presented in this 
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Fig. 2. — The dipole {1 = 1) normalized growth rates <n (black 
dots connected with continuous lines) in terms of the pulsation pe- 
riods for a 0.530A/q PG1159 template model located before the 
evolutionary knee (T eff = 138400 K, log(L*/L Q ) = 3.14). Num- 
bers indicate the radial order k for low order modes. The large 
numerical range spanned by r] is appropriately scaled for a better 
graphical representation. The two ranges of overstable g-modes 
— one due to the K-mechanism and the other induced by the e- 
mechanism — are clearly discernible. Plus symbols connected with 
dotted lines correspond to the case where the e-mechanism is ex- 
plicitly suppressed in the stability calculations. 

work was carried out with the linear, nonradial, nonadi- 
abatic pulsation code described in Corsico et al. (2006). 
The "frozen-in convection" approximation was assumed 
because the flux of heat carried by convection is neg- 
ligible in PG1159 stars. At variance with Corsico et 
al. (2006), in this work we have fully taken into ac- 
count the e-mechanism for mode driving operating in the 
He-shell nuclear-burning region. Because we are inter- 
ested in PG1159 stars which are H-deficient, we are only 
concerned with the He-burning reactions. Fortunately, 
because g-mode pulsation timescales are much shorter 
than the timescales of nucleosynthesis, possible phase 
delays between the temperature perturbations and the 
abundance variations are unimportant. Hence, they can 
be neglected, largely simplifying the pulsational stability 
analysis (Unno et al. 1989; Kawaler et al. 1986) 2 . 

3. NONADIABATIC RESULTS 

We analyzed the stability properties of about 4000 
stellar models covering a wide range of effective tem- 
peratures (5.5 > log(Tcff) > 4.7) and stellar masses 
(0.515 < M*/M < 0.741). For each model we restricted 
our study to £ = 1 g-modes with periods in the range 
50- 7 000 s. 

2 We note that some of our sequences have trace surface abun- 
dances of H (Xn < 10 — 3 ) which give rise to some H burning. 
However, exhaustive test stability computations demonstrate that 
H burning is very weak and extends on a extremely narrow layer, 
as a result of which the H-shell burning is completely irrelevant 
in destabilizing modes and will not be further considered in this 
paper. 



3.1. A single template model 

We start our description by focusing on a 0.530M© 
PG1159 template model with T off = 138 400 K and 
log(i»/L Q ) = 3.14 located before the evolutionary knee 
in Fig. [T] The surface chemical composition of the 
model is X( 4 He) = 0.33, X( 12 C) = 0.39, X( 13 C) = 0.05, 
X( U N) = 0.02, and X{ 16 0) = 0.17. Fig. [2displays the 
normalized 1=1 growth rates r\ — — 3?(<t)/Sj((t) (where 
5ft(er) and 3(c) are the real and the imaginary parts, re- 
spectively, of the complex eigenfrequency a) in terms of 
the pulsation periods (n) corresponding to our template 
model. In the interests of a better graphical representa- 
tion, the huge numerical range spanned by r\ is appropri- 
ately scaled (see Gautschy 1997). The sign function al- 
lows to discriminate between stable and unstable modes. 
The presence of two well-defined families of overstable g- 
modes, one at the intermediate- and long-period regime, 
and the other one at the short period regime, is apparent. 
The first group of periods (« 750 — 1800 s) corresponds 
to modes driven by the well-known K-mechanism oper- 
ating at the region of the opacity bump due to partial 
ionization of C and O, centered at logT w 6.2 (Gautschy 
et al. 2005, Corsico et al. 2006). The second group of pe- 
riods, which are associated to low radial order g-modes, 
are destabilized by the action of the vigorous He-shell 
burning through the e-mechanism. 

The short-period instabilities uncovered here are of the 
same nature than those found by Kawaler et al. (1986) in 
the context of H-deficient hot central stars of planetary 
nebulae. Here, as in that work, the e-mechanism induced 
by the He-shell burning constitutes the source of driv- 
ing. In absence of this destabilizing agent, the overstable 
modes with periods in the range (« 100— 180 s) turn out 
to be stable, while the remainder modes of the pulsation 
spectrum remains unchanged. This is vividly displayed 
in Fig. [21 that shows with plus symbols the results of ad- 
ditional stability computations in which the nuclear en- 
ergy production rate, e, and the logarithmic derivatives 
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in the pulsation equations. It is worth emphasizing that 
in the present effort we are able to obtain destabilization 
of g-modes through both the K-mechanism and the e- 
mechanism in the same PG1159 equilibrium model. This 
is at variance with the study by Kawaler et al. (1986), 
who reported only e-destabilized modes. 

The e-mechanism behaves as an efficient filter of 
modes that destabilizes only those (/-modes that have 
their largest maximum of the temperature perturbation 
(ST/T) in the narrow region of the He-burning shell (see 
Kawaler et al. 1986). In the specific case of our template 
model, only the g-modes with k = 3,4,5, and 6 meet 
such a condition and, as a result, they are e-destabilized. 
For k = 1, 2 the largest maximum of ST/T lies at inner 
layers with respect to the He-burning shell. Thus, these 
modes are stable. For modes with k > 7 the opposite is 
true and these modes also are stable. 

Test stability calculations with 1 = 1 for our template 
model indicate that there exist only one quadrupole e- 
dcstabilizcd g-mode which corresponds to k — 5 with 
a period n ~ 95 s, about 40% shorter than the corre- 
sponding k — 5 mode with I = 1 (n ~ 157 s). Hence, in 
general, a narrower range of shorter periods is expected 
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to be associated with e-destabilized (/-modes with £ = 2 
as compared with the case of I = 1 . 

3.2. A new PG1159 instability strip 

Having described our results for a single template 
model, we now are in conditions to examine the lo- 
cation and extension of the complete unstable domain 
associated with the e-mechanism. Our results are de- 
picted in Fig. [3J which displays the instability strip 
of e-destabilized modes in the \ogT e Q — log g drawn 
with thick black curves along the PG1159 evolutionary 
tracks. In addition, the GW Vir instability domain of re- 
destabilized modes (see Corsico et al. 2006) is depicted 
with red (gray) lines along the tracks. Note that the 
instability strip for e-destabilized modes partially over- 
laps the domain of K-destabilized modes. So, our results 
indicate the existence of three well-defined instability 
regimes in the logT c ff — \ogg plane: a regime — splittcd 
into two regions, one at low gravity and the other at 
high gravity — in which stellar models harbour interme- 
diate/long period g-modes excited by the k- mechanism 
only, another one corresponding to short-period modes 
destabilized by the e-mechanism only, and finally a re- 
gion in which models experience pulsational destabiliza- 
tion by the K-mechanism and the e-mechanism of driving 
simultaneously Notably, the region corresponding to the 
e-mechanism only is not occupied by any known PG1159 
star (see Fig. [3]). 

We stress that in previous works (Kawaler et al. 1986; 
Saio 1996; Gautschy 1997) only some short-period g- 
modes were found to be destabilized by the e-mechanism. 
Needles to say, due to the very few e-destabilized modes 
found in those exploratory works, no clear extension and 
location of the e-mechanism instability domain were ob- 
tained, thus hampering those authors from making fur- 
ther consideration of such modes. At variance with those 
works, in the present study we are able to find a complete 
instability strip of e-destabilized modes. 

The degree of driving, and the place that it might oc- 
cur in the logT e ff — log g diagram, is sensitive to the 
stellar mass, previous evolutionary history, and so on. 
Thus, due to the uncertainties in the stellar evolution 
modelling (overshooting, nuclear reaction rates, etc), the 
surface and internal composition of PG1159 stars are not 
known in detail, and so a clear instability domain for e- 
dcstabilizcd pulsations is difficult to drawn. So, what 
is shown in the Fig. [3J is the shape and location of the 
e-mechanism instability strip obtained by us under the 
particular assumptions adopted in the construction of the 
PG1159 evolutionary models of Miller Bertolami & Al- 
thaus (2006). The extension and location of this instabil- 
ity domain might change if other assumptions for the evo- 
lutionary history of the progenitor stars were adopted. 

All of the overstable e-destabilized modes computed in 
this work are characterized by very tiny (10~ 9 — 10 -12 ) 
linear growth rates n, by far smaller than those char- 
acterizing overstable modes excited by the K-mechanism 
(10 -8 < n < 10~ 4 ). So, the question rises about what 
would be the chance for a given e-destabilized mode to 
have plenty of time for developing observable amplitudes. 
To analyze this question we consider the time interval 
that the models spend crossing the instability strip of 
e-destabilizcd modes, At, and the maximum and mini- 
mum e-folding times r° lax and r™ m , respectively, of the 




Fig. 3. — Same as Fig. [T] but including the loci of models 
having 1=1 (dipole) re-destabilized modes with solid red (gray) 
curves along the tracks, and models harbouring short-period e- 
dcstabilized modes according to the present study. Superposition 
of both types of curves corresponds to stellar models with both e- 
and re-destabilized modes (shaded area). The location and desig- 
nation of relevant PG1159 stars is also shown [Color figure only 
available in the electronic version of the article] . 

unstable modes for a given stellar mass. The e-folding 
times are defined as r e = l/|5(er)|, such that the time 
dependence of the amplitude of the pulsations is given 
by £(t) oc e l(Tt , and $s(<r) < for overstable modes. 

The values of At, T™ ln , and r™ ax are provided in Table 
[T]for each value of the stellar mass. Note that the three 
timescales monotonically decrease for increasing stellar 
mass. For all of our PG1159 sequences we found that 
the most unstable modes — those with the smaller values 
of r e — are found near the low-gravity (high-luminosity) 
boundary of the instability domain (upper black dashed 
line in Fig. [3]), when the models are still evolving to 
the blue before reaching the evolutionary knee. On the 
contrary, when models are already evolving toward the 
white dwarf cooling track, the e-destabilized modes are 
only marginally unstable, and so they are characterized 
by large e-folding times. 

Table [T] shows that r™ ln <C At for all of our sequences. 
This means that g-modes that are destabilized at epochs 
before the evolutionary knee, characterized by short e- 
folding times, have time enough to reach observable am- 
plitudes before the star leaves the instability strip. On 
the other hand, it is apparent that r™ ax > At. Thus, 
the g-modes that are destabilized in models close to the 
high-gravity limit (low-luminosity) of the instability strip 
(lower black dashed line in Fig. [3]) have little — or even 
null — chances to develop observable amplitude before 
the model abandons the instability domain. 

In summary, our computations predict that some g- 
modes (those with short r e ) could have plenty of time to 
grow and finally develop observable amplitudes. We cau- 
tion, however, that this prediction is based on a linear 
stability analysis, and that the last word should came 
from a detailed non-linear description of nonadiabatic 
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TABLE 1 

THE MINIMUM AND MAXIMUM e-FOLDING times (in yr), and the 
TIME (IN YR) THAT PG1159 MODELS SPEND WITHIN THE 
INSTABILITY STRIP OF e-DESTABILIZED MODES. 



M./M© 


min 
' e 


max 
' e 


At 


0.515 


3410 


1.5 x 10 6 


1.60 x 10 5 


0.530 


2580 


1.0 x 10 6 


1.01 x 10 5 


0.542 


1610 


3.8 x 10 5 


5.95 x 10 4 


0.565 


1400 


1.3 x 10 5 


2.78 x 10 4 


0.589 


1160 


1.0 x 10 5 


2.47 x 10 4 


0.609 


742 


4.7 x 10 4 


1.26 x 10 4 


0.664 


361 


1.8 x 10 4 


4830 


0.741 


180 


7000 


1570 



pulsations. Such a nonlinear treatment is not available 
at the present stage. Also, there are other effects (stel- 
lar winds, diffusion, etc) suspected to be present in real 
stars, that could be affecting the growth of pulsations. 
The assessment of their effects on the modes predicted 
to be unstable in the frame of our analysis is beyond the 
scope of the present study. 

The next step in our analysis is to derive the range of 
periods (II) of overstable e-destabilized modes. Fig. 2] 
displays the regions of the K-mechanism instability do- 
mains in the log Toff — II diagram, depicted with lines 
of different colours for the various stellar masses. No- 
tably, the figure also shows the presence of a separate, 
well-defined instability domain for a broad range of ef- 
fective temperatures (5.46 > logT e g > 4.99) and pul- 
sation periods in the interval 55 < II < 200 s, associ- 
ated to e-destabilized g-modes with radial orders rang- 
ing from 2 to 5 for M* = 0.5I5M Q and from 3 to 8 for 
M* = O.741M . The stages corresponding to phases be- 
fore (after) the evolutionary knee are depicted with small 
dot (plus) symbols. A close inspection of the figure re- 
veals that for the low-mass models, most of modes are 
destabilized after the evolutionary knee. For the high- 
mass models the situation is reversed, that is, most of 
overstable modes are destabilized before the evolution- 
ary knee. The existence of this new instability domain 
of short-period g-modes in stellar models representative 
of PG1 159 stars is the main result of our study. 

In particular, it is worth emphasizing that the e- 
mechanism should be active in a PGII59 star irrespective 
of the precise abundances of He, C, and O at the surface, 
because in this case the mode excitation takes place at 
deep layers in the star. This is at variance with the k- 
mechanism, which is strongly dependent on the exact 
O/C/He abundances at the driving regions (see Quirion 
et al. 2007). 

3.3. The case of the PNNV star VV 47 

An immediate prediction of the present study is that 
PGII59-type stars populating the overlapping region of 
Kr and e-destabilized modes in the log T c g — log g diagram 
(the shaded region in Fig. [3]) should exhibit both short- 
and intermediate/long-period luminosity variations si- 
multaneously. Tablc[5]lists the known PGII59 candidate 
stars. A glance of this table leads us to a somewhat dis- 
appointing conclusion: most of the stars located in the 
region of interest are not variables at all or have not been 
scrutinized for variability. Other stars, at most, exhibit 
intermediatc/long-period luminosity variations which are 
typical of the high/intermediate-order g-modes driven 
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Fig. 4. — The dipole [£ = 1) instability domains for overstable 
K-destabilized g-modes, shown with thin lines of different colours 
for the various stellar masses. The e-mechanism instability domain 
is emphasized with a shaded area. Short-period dipole unstable 
e-destabilized g-modes are depicted with dot (plus) symbols for 
stages before (after) the evolutionary knee. Numbers indicate the 
radial order of the modes. [Color figure only available in the elec- 
tronic version of the article] . 

by the K-mechanism, but not the expected short peri- 
ods typical of e-destabilized modes. In particular, this is 
the case for the prototype DOV star, PG 1159—035. 

There is one object, the PNNV star VV 47, which is 
suspected to pulsate in short- and long-period modes 
(Gonzalez Perez et al. 2006). This star (T off = 130 000 ± 
5000 K, log,? = 7 ±0.5, C/He= 1.5 and 0/He= 0.4; 
Werner & Herwig 2006) was first observed as potentially 
variable by Liebert et al. (1988). Later, it was monitored 
by Ciardullo & Bond (1996), but no clear variability was 
found. Finally, Gonzalez Perez et al. (2006) were able 
to confirm the — until then, elusive — intrinsic variabil- 
ity of VV 47 for the first time. They found evidence that 
the pulsation spectrum of this star is extremely complex. 
The most outstanding feature of VV 47 is the presence 
of high-frequency peaks (at periods ~ 130 — 300 s) in the 
power spectrum, which could be serious candidates for 
low-fc radial order g-modes triggered by the e-mechanism. 

We decided to test the attractive possibility that the 
short-periods observed in VV 47 could be due to the 
e-mechanism. We first estimated the stellar mass of 
VV 47. From the location of VV 47 in the log T e g — log g 
plane (see Fig. [T|) it is apparent that the spectroscopic 
mass of VV 47 is of re 0.525A/©. In addition, a pre- 
liminary adiabatic asteroseismological analysis on this 
star indicates that the seismological mass of VV 47 - 
obtained from the period spacing data (All re 24 s) — 
is of re 0.52 — 0.53A1©, in excellent agreement with the 
spectroscopic derivation. So, we shall focus on the case 
of the evolutionary sequence of M* = 0.530A1©. This 
sequence is characterized by a thick He-rich envelope 
(M env ~ 0.045M©). 

We would like to see how well the theoretical ranges 
of periods of unstable modes corresponding to this se- 
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TABLE 2 

Known PG1159 stars populating the overlapping 

INSTABILITY REGION OF €- AND re- DESTABILIZED MODES. 



Period range [s] 





Star 


PN 


Variable 


1 


FG 1159—035 


no 


yes 


2 


NGC 650—1 


yes 


no 


3 


VV 47 


yes 


yes (.') 


4 


PG 1144+005 


no 


no 


5 


Jn 1 


yes 


yes (?) 


6 


Abell 21 


yes 


no 


7 


K 1-16 


yes 


yes 


8 


Longmore 3 


yes 


no 


9 


HS 2324+3944 


no 


yes 


10 


HS 1517+7403 


no 


no 


11 


PG 1424+535 


no 


no 


12 


IW 1 


yes 


no 


13 


Sh 2-68 


yes 


? 


11 


HS 0704+6153 


no 


? 



430 - 840 

~ 260 
~ 2170-4300 

454 - 1860 

1500 - 1700 

2005 - 2570 



quence fit the observed period spectrum of VV 47. Fig. 
[5] displays the regions of the K-mechanism instability do- 
main (light and dark grey) for the 0.530M Q sequence. 
The figure also shows the presence of a well-defined in- 
stability domain (77 < II < 180 s) that corresponds 
to e-destabilized <?-modes with k = 2, • • • , 6 (large and 
small dots). Also depicted in the plot are the periods 
reported by Gonzalez Perez et al. (2006) for VV 47 with 
their corresponding uncertainties in T e ff. We have em- 
phasized with black small circles the periods associated 
with modes having the best chances to be real, accord- 
ing to Gonzalez Perez et al. (2006). It is apparent that, 
whereas most of the long periods observed in VV 47 are 
qualitatively explained by the K-mechanism when the 
model star is before the evolutionary knee, the short- 
period branch (below ~ 300 s) of the pulsation spectrum 
of the star is not accounted for at all by the theoretical 
domains corresponding to this destabilizing agent. We 
can see, instead, that the short periods of VV 47 — in 
particular II < 200 s — are satisfactorily accounted for 
by the e-destabilized g-modes. Note, however, that if 
only periods detected with sufficiently high significance 
(black filled circles) are used to compare with our the- 
oretical predictions, then the period at 261.4 s cannot 
be explained by a low-order <?-mode excited by the e- 
mechanism. In fact, this period is considerably longer 
than the longest period (« 180 s) of the g-modes which 
can be excited by the e-mechanism as our analysis pre- 
dicts. 

Clearly, more observational work is needed to put the 
reality of the short periods in VV 47 on a solid basis. 
Were the existence of these short periods confirmed by 
future observations, then they could be attributed to the 
e-mechanism, and this could be indicating that W 47 
should have a thick He-rich envelope as to support an 
active He-burning shell. 

3.4. The case of the DOV star PG 1159-035 

Another consequence of our investigation concerns the 
pulsating star PG 1159—035, the prototype of the class 
and the best-studied DOV. Indeed, note from Fig. Othat 
a trend of our results is that this variable star should ex- 
hibit short-period e-destabilized modes if the thick He- 
rich envelopes derived from our evolutionary calculations 
were representative of the star. These modes are not ob- 
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Fig. 5. — The I = 1 regions of the K-mechanism instability do- 
main, shown with light (dark) grey for stages before (after) the 
evolutionary knee, corresponding to the 0.530AfQ sequence. Ar- 
rows indicate the time sense of evolution. Also shown is the evo- 
lution of the periods corresponding to the e-destabilized modes 
92, • • • ,96: with large (small) dots for stages before (after) the evo- 
lutionary knee. Finally, the periods reported by Gonzalez Perez et 
al. (2006) for VV 47 with their corresponding uncertainties in T e ff 
arc displayed with small circles. Periods detected with sufficiently 
high significance are emphasized with black filled circles. 

served by Costa et al. (2007). This result suggests that 
the He-burning shell may not be active in PG 1159—035. 
This would indicate that this star has a thinner He- 
rich envelope than what is traditionally derived from 
standard evolutionary calculations, in line with the re- 
cent finding by Althaus et al. (2008) that a thinner 
He-rich envelope (at least a factor of two below of the 
value predicted by the standard evolution theory) for PG 
1159—035 should be invoked to alleviate the longstand- 
ing discrepancy between the observed (Costa & Kepler 
2008) and the theoretical (Corsico et al. 2008) rates of 
period change in that star. 

If the short periods observed in VV 47 were confirmed, 
then we should face the problem of the coexistence of 
two PG1159 stars located very close each other in the 
logT c ff — logg diagram (see Fig. [3]) but with substan- 
tially different thickness of the He-rich envelopes. This 
would suggest that these stars could have had a differ- 
ent evolutionary history, a suggestion reinforced by the 
fact that VV 47 still retains a planetary nebula while PG 
1159-035 does not. 

4. SUMMARY AND CONCLUSIONS 

In this paper, we have presented a fully nonadiabatic 
stability analysis on state-of-the-art PG1159 models gen- 
erated taking into account the complete evolution of pro- 
genitor stars, through the thermally pulsing AGB phase 
and born-again episode. We have explored the possi- 
bility that nonradial <?-mode pulsations could be desta- 
bilized by a He-burning shell through the e-mechanism. 
Our study covers a broad range of stellar masses and 
effective temperatures. We confirm and extend the pio- 
neering work of Kawaler et al. (1986), Saio (1996) and 
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Gautschy (1997) on this topic. The main results are the 
following: 

- We found strong theoretical evidence for the exis- 
tence of a separate, well-defined PG1159 instability 
strip in the log T e g — log g diagram characterized by 
short-period g-modes excited by the e-mechanism 
due to the presence of active He-burning shells. No- 
tably, this instability strip partially overlaps the al- 
ready known GW Vir instability strip due to the 
K-mechanism acting on the partial ionization of C 
and/or O in the envelope of the PG1159 stars. We 
emphasize that while in previous works only some 
short-period g-modes were found to be destabilized 
by the e-mechanism, in the present study we found 
a complete instability strip of e-destabilized modes. 

- At variance with the classical /c-mechanism respon- 
sible for the intermediate/long-period GW Vir pul- 
sations, the e-mechanism should be efficient even 
in PG1159 stars with low C and O content in their 
envelopes. 

- The e-driven g-modes that are destabilized at 
epochs before the evolutionary knee are character- 
ized by short e-folding times (between « 180 yr for 
M* = 0.741M Q and w 3000 yr for M* = O.515M ), 
and so they probably have time enough as to reach 
observable amplitudes before the star leaves the in- 
stability strip. Note, however, that nonlinear ef- 
fects, or the presence of a variety of phenomena 
such as stellar winds or diffusion, could affect the 
growth of pulsations. 

- We have closely examined the case of VV 47, 
the only PG1159 star for which observational evi- 
dence of the presence of short-period g-modes ex- 
ists (Gonzalez Perez et al. 2006). For this star we 
have derived for the first time a seismological mass 
of « 0.52— 0.53Mq, in excellent agreement with the 
spectroscopic mass (« O.525M ). If we accept that 
all of the periods reported by Gonzalez Perez et al. 
(2006) are real, our stability analysis provides very 
strong support to the idea that the physical origin 
of the short periodicities could be the e-mechanism 
powered by an active He-burning shell, whereas the 
long-period branch of the period spectrum of this 
star should be due to the K-mechanism acting on 
the region of partial ionization of C and O. How- 
ever, if only periods detected with sufficiently high 
significance are taken into account, then the pe- 
riod at 261.4 s can not be explained by a low-order 
g-mode excited by the e-mechanism. 



We speculate that the absence of short periods (< 
300 s) in the pulsation spectrum of PG 1159—035 
could be indicating that the He-burning shell may 
not be active in this star. This would indicate that 
PG 1159—035 has a thinner He-rich envelope than 
what is traditionally derived from standard evolu- 
tionary computations. 



In light of our results, if the reality of the short peri- 
ods of VV 47 were confirmed by follow-up observations, 
this star could be the first known pulsating PG1159 
star undergoing nonradial g-modes destabilized by the e- 
mechanism. Even more, VV 47 could be the first known 
pulsating star in which both the ft-mechanism and the 
e-mechanism of mode driving are simultaneously oper- 
ating. Further time-series photometry of VV 47 will be 
needed to firmly establish the reality of the short-period 
pulsations detected in this star. 

On the other hand, the apparent absence of short- 
period pulsations in the remainder variable stars — such 
as K 1-16, HS 2324+3944, and Jn 1— could be an indi- 
cation that, like PG 1159—035, they are characterized by 
thin He-rich envelopes, as a result of which they should 
lack of stable He-shell burning. Another possibility is 
that short-period pulsations could be indeed present in 
these stars, but with very low amplitudes, below the ac- 
tual detection limits. 

Also, it is quite intriguing the absence of both short- 
and intermcdiate/long-pcriod pulsations in the up to 
now constant stars (NGC 650-1, PG 1144+005, Abell 
21, Longmore 3, HS 1517+7403, PG 1424+535, IW 1) 
that populate the overlapping region of the e- and k- 
dcstabilizcd modes. In any case, extensive searches for 
low amplitude intrinsic variability in these stars and also 
in the stars Sh 2-68 and HS 0704+6153, which have not 
been observed for variability yet, should be worth do- 
ing in order to test the existence of the new theoretical 
instability strip uncovered in this work. 
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